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The reaction pp — > ppTr" was studied with the WASA detector at the CELSIUS storage ring. The 
center of mass angular distribution of the tt" was obtained by detection of the 7 decay products 
together with the two outgoing protons, and found to be anisotropic with a negative second derivative 
slope, in agreement with the theoretical predictions from a microscopic calculation. 
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I. INTRODUCTION 

The first high precision measurements of single neutral 
pion production in nucleon-nucleon collisions, using stor- 
age ring technology, were done more than a decade ago. 
Still, the theoretical interpretation of the dominant pro- 
duction mechanism remains uncertain. The magnitude of 
the total cross section of the pp ppn^ reaction in the 
threshold region, where only angular momenta equal to 
zero arc important in the final state, was measured to 
be about five times larger 2, 3] than what was predicted 
by the theoretical models available at the time. However, 
the energy dependence was found to be consistent with 
the widely accepted Koltun and Reitan model |j| based 
on s-wave pion production and rescattering The 
experimental result was confirmed and expanded even 
closer to threshold Q whereas the large theoretical ac- 
tivity that was triggered by the new high precision data, 
brought conflicting not yet settled results. 

The first successful remedy to fill in the discrepancy 
between experiment and theory was to take into account 
the exchange of heavy mesons Q . The off-shell pion 
rescattering (together with the Born term) was also sug- 
gested to fill in the gap in the cross section Both 
these theories cannot be right unless there are some other 
additional effects. Furthermore, approaches using chiral 
perturbation effective field theory (ChPT) reached a dif- 
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ferent conclusion than meson field theory and found the 
interference between the direct term and the pion rescat- 
tering to be destructive 0, 0. Improved calculations 
carried out in momentum-space increased the rescatter- 
ing arnplitude for the ChPT treatment by a factor of 
three [ij. Considerable progress has since been made 
developing the ideas of , using an ordering scheme that 
takes into account the large momentum transfer typical 
for meson production in NN collisions [13 . [T3. [l4 . J^J . 
Within this scheme it was possible to describe the reac- 
tion pp — > dTT"*" near threshold [T^ and a corresponding 
study of the 7r° production is under way [T^ . 

A calculation taking into account the exchange of 
two different heavy mesons, pion rescattering and the 
Pii(1440) nucleon resonance reproduced the total cross 
section numbers [l^- Relativistic effects were studied in 
the impulse approximation in Ref. |19|. The exchange 
of the mesons tt, p, uj and a, with the nucleon and the 
A(1232) isobar as intermediate states, using a relativis- 
tic treatment in a covariant one-boson exchange model 
over an energy range from near threshold to 2 GeV, gave 
reasonable agreement with data The effect of the 

resonances Pii(1440), 511(1535) and L»i3(1520) together 
with the impulse and pair diagram terms were studied in 
Ref.Hll 

The possible influence on the differential cross sections 
due to contributions from higher partial waves, inclu ding 
(i- waves, was studied experimentally at CELSIUS |23| . 
Angular distributions as well as total cross sections were 
recently measured from threshold up to 10 MeV above, 
by the TOF collaboration using an extracted beam [23l |. 
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FIG. 1: Left panel : Scatter plot of the two proton polar an- 
gles in the laboratory system, obtained from a monte carlo 
simulation of phase space distributed pp — > ppix^ events. The 
geometrical coverage of the WASA detector is shown by the 
overlaid histogram. Right panel : The relative proton mo- 
menta {p from a simulation according phase space. The la- 
belled curves represent different ranges of p corresponding 
to their detection in a full-blown simulation of the detector 
setup: two forward prongs (2FD), one forward and one central 
(IFD-flCD). The two trigger conditions used in the experi- 
ment are of the types 2FD and IFD + ICD (see sec. Ill At . 



The angular distributions were isotropic as expected close 
to threshold. The total cross section obtained was about 
50 % larger than the lUCF 1] and CELSIUS data [|. 
The reason for the large deviation compared to the pre- 
vious storage ring experiments was suggested to be due 
to a significant loss of events in the internal experiments, 
where the very forward going protons escape down the 
beampipe undetected. At threshold a strong final state 
interaction could cause the loss of a large number of pro- 
tons that would not be properly accounted for. 

The data set was drastically increased when the reac- 
tion — > ppTT^ was measured at beam energies be- 
tween 325 and 400 MeV by the PINTEX collaboration 
at lUCF ^ 25, 26, 27]. All possible polarization ob- 
servables were deduced in the kinematically complete ex- 
periment and a general formalism was developed from a 
partial wave analysis, in order to obtain an expansion of 
the observables in terms of a complete set of functions 
mapping the angular dependence |23. Thus an analy- 
sis method was realized with all the physics information 
contained in the deduced coefficients. However, contri- 
butions from Ds and Sd to the final state were not taken 
into account. The only theoretical model that so far has 
been compared to these data is the microscopic model de- 
veloped by the Jiilich group ,2&j_29j. The phenomenology 
of the model includes direct production, s- and p-wave 
rescattering of the pion, pair diagrams and excitation of 
the A(1232). Angular momenta up to Lp,lq < 2 be- 
tween the two protons and the pion with respect to the 
NN subsystem, are included. The same group has re- 



cently performed a partial wave analysis using the data 
and the assumptions of j23| and compared the extracted 
quantities to those of their meson-exchange model [30l| . 
Most of the amplitudes are shown to be reproduced fairly 
well by the model, except for the amplitude ^Pi '^PqP 
that deviates significantly from what is extracted from 
the data. For a quantitative assessment of pp — > ppir'^ it 
also turns out that the A excitation plays a major role. 
For a summar y on near threshold meson production ex- 
periments s ee ^11 . The status of the theoretical field is 
reviewed in |32j . 

In spite of the interest in the reaction pp — > ppn'^ dur- 
ing the last 15 years, the reports on the pion angular 
distributions at energies up to 400 MeV, suffer from low 
statistics and/or small acceptance. We have measured 
the unpolarized angular cross section of the 7r° at 400 
MeV, with the aim to resolve some of the ambiguities 
and complement the partial wave analysis based solely 
on polarized data. 



II. MEASUREMENT 

The experiment was done using the WASA 47r detec- 
tor facility ^] situated in the CELSIUS accelerator and 
storage ring at Uppsala, Sweden. A stored circulating 
proton beam of energy 400 MeV was let to interact with 
a stream of small {(f) ^ 30 fim) frozen hydrogen pellets. 
All three outgoing particles from the reaction pp — > ppir'^ 
were detected. The protons were fully stopped either 
both in the Forward Detector (FD), or one in the FD 
and one in the Central Detector (CD). The FD consists 
of a stack of scintillator and wire chamber planes, pri- 
marily adapted to measure the four momenta of recoil- 
ing nuclei. The CD is constructed for measuring meson 
decay products, and comprises the Scintillating Electro- 
magnetic Calorimeter (SEC) made up of 1012 Csl detec- 
tor elmeents, a Plastic Barrel (PB) for charged particle 
detection and a Mini Drift Chamber (MDC) for measur- 
ing the momenta of charged particles. In the current ex- 
periment only the FD and the SEC were used for energy 
measurement and the PB for the rejection of charged 
particles. Since the energy and angular resolution for 
protons was relatively poor in the CD, the final analysis 
was based on the detection of the tt'^ 27 decay in the 
CD. 



A. Data Analysis 

The event selection was handled using two different 
sets of criteria based on two different track types, which 
were either both protons detected in the FD {2FD-type), 
or one in the FD and one in the CD {IFD ICD -type). The 
requirements were coincident fast signals from either two 
hits in one scintillator layer of the FD or one in the FD 
and one hit in the forward part of the PB. These triggers 
gave an unbiased acceptance of the CD but yielded very 
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FIG. 2: The invariant mass of of the two 7's selected by 
the IFDl CD-type trigger. Left panel : The cut based on 
the constraint that the sum of the energies of the two 7's is 
within the kinematical limits for vr" production, is shown by 
the dashed line. The bold line depicts what is cut out based 
on the relation between the opening angle and the planarity 
of the two 7's, representing largely background from elastic 
scattering. Right panel : The combination of the two cuts 
is shown by the dashed line and the final resulting invariant 
mass of the two 7's, is drawn by the bold line. 
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FIG. 3: The experimental center-of-mass tt angular distribu- 
tion, arbitrarily normalized. The shaded area is the result of a 
simulated phase space generated isotropic distribution of the 
pp ppiv'^ after passing through the detector system. The 
solid line corresponds to the predicted histogrammed values 
from a simulation weig hted with the microscopic calculation 
by Hanhart et al. |28L l29|. The statistical uncertainties are 
negligible. 

high countrates, why prescaling was necessary. The main 
background reactions present in the trigger were pp — > 
dTT"*" and pp elastic scattering. 

The geometrical acceptance for detection of the out- 
going protons from the reaction pp — > ppn'^ is shown in 
Fig.n The angular coverage was 3° - 17° and 20° - 155°, 
for the FD and the CD respectively. Since there were no 
triggers set for the case when both protons are emitted 
at 9iab > 17°, these events escape the current analysis. 



However, the full range of the relative momenta of the 
two protons is covered by the experiment, (c.f. the right 
panel of Fig. ^| , which is crucial from the physics inter- 
pretation point-of-view. 

The basic condition for an accepted event of the 
2FD — type was particle identification of the protons in 
the FD done by AE — E technique, and the presence in 
the CD of two neutral tracks from the decay of tt" — > 27. 
Additional constraints were based on the comparison of 
the reconstructed polar and azimuthal laboratory angles, 
plus cuts in the center-of-mass energy, with respect to the 
missing mass of the two protons and the invariant mass 
of the two 7's from the 7r° decay. The conditions ap- 
plied were: \9mx - Oim\ < 15°, \(j)Mx ~ < 15° and 
\Emx — Eim\ < 30 MeV, where Mx is the missing mass 
of the two protons and IM the invariant mass of the two 
7's. The consistency of the tt*^ angle reconstructed from 
the missing mass and the invariant mass, respectively, 
was investigated. The two approaches agreed and thus 
the conclusion was drawn that analysis of the IFDl CD- 
type prongs could be done using CD information only. 

The selection of event candidates for the type of events 
with one forward and one central prong {IFDICD — 
type) was done by particle identification of the FD pro- 
ton. Furthermore, for the two 7's a cut was based on 
the relation between the opening and planarity angles. 
An additional constraint applied was that the sum of the 
energies of the two 7's was within the kinematical lim- 
its for 7r° production, i.e. 135 < EE^ < 238 MeV, The 
TT*^ peak obtained from the invariant mass of the two 7's, 
before and after track requirements are fulfilled, is shown 
in Fig.ia 

The two sets of selected ppir^ events were weighted to- 
gether according to their relative trigger prescaling fac- 
tor. The experimental angular distribution of the tt'' in 
the center-of-mass, uncorrected for acceptance, is seen 
in Fig. 13 Displayed are also simulations using either 
isotropically distributed events according to phase space 
or events weighted with the theoretical calculation of the 
Jiilich model by Hanhart et al. |2^ I^Pl . More details on 
the data reduction procedure can be found in |3^ . 



III. RESULTS 

The acceptance corrected 7r° angular distribution is 
shown together with the prediction by the Jiilich model 
in Fig. 1^ The experimental data points and the theoret- 
ical curve are normalized to atot — 93 ± 7.2 /ibarn from 
|22j| . The systematic uncertainties dominate, primarily 
emanating from the acceptance varying with the central 
detector's geometry. In order to obtain an estimate of 
the magnitude of this error the outmost layers in the for- 
ward and the backward parts respectively, were excluded 
in the analysis. 

In previous experimental reports p3 . l35ll3^l3^l33 . l3^ 
liol IIH concerning the angular distribution a slope pa- 
rameter b was defined according ^ ^ -\- b cos^ 67^ , see 



4 



10 



u 

I 

OS 

a 

e5 





full range of p 


- 














• Wasa 




• Hanhart 



0.5 



-1 





COS (6 sC) 



FIG. 4; Acceptance corrected center-of-mass tt" angular dis- 
tribution, normalized to otot ~ 92.3 ± 7.2 /ibarn |2^. The 
uncertainites are dominated by systematic effects. The solid 
line corresponds to the microscopic calculation by Hanhart et 
al. |2^ . |29( . The dashed line represents a fit of the dependence 
on cos'^ d-n, taking only the statistical errors into account. 
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FIG. 6: Acceptance corrected center-of-mass tt" angular dis- 
tribution, for p < 53 MeV/c, compared and normalized to 
[2^ . The solid line represents a fit of the WASA data points 
to the dependence on cos^ 6,r, taking only the statistical errors 
into account. 



TABLE I; Correspondence between Iri -|- 7 and b 
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FIG. 5: A compilation of measurements of the slope param- 
eter below 1 GeV is shown. The definition of 6 differs by a 
factor of three in |25| why the values given by that reference 
have been divided by three for consistency. (The large error 
bar at 400 MeV is from )35|l ). 



FigEl for a compilation of measurements below 1 GeV. 
There is a large spread of the values, probably mainly 
due to the varying cov erag e of p in the measurements. 
One recent experiment |22l | yield a negative h up to 360 
MeV beam energy but at 400 MeV the slope was re- 
ported to be positive in discrepancy with the present re- 
sult h — —0.116 ± 0.010. However the acceptance of [S^l 
was limited with respect to p, with a coverage similar 
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to the 1FD — type case, and a model dependence was 
introduced by extrapolating into unmeasured regions of 
phase space. It should be noted that events with both 
proton angles larger than 17° were not detected within 
the current acceptance, see Fig. ^ There are indications 
that the inclusion of such events would slightly flatten 
the distribution 42] . Whereas the 7r° angular distribu- 
tion integrated over all -p displays discrepancies among 
the different experiments, a selection of S'-wave protons 
might shed some light. At very low relative momenta 
between the protons, p < 53 MeV/c, the two CELSIUS 
measurements agree at least qualitatively, see Fig. El At 
800 MeV beam energy using the same cut in o, an even 
larger negative second derivative was found 43], which 
was also predicted by a phenomenological calculation 44] . 

A direct comparison between the present experiment 
and the expansion deduced from the double polarization 
data |27jl can be done by using the sampling method ,4j] , 
that allows to integrate the prediction of a theoretical 
model over the experimentally accessible phase space re- 
gion, see Fig. [3 The coefficient combination + / 
and the dependency on the cos^(07r)-term (Eq. 11 27]), 
have been taken into account, all other variables are ig- 
nored. See Table for the correspondence between the 
slope parameter b and H^'^ + 1. 

The prospect of improving the accuracy of certain co- 



5 




-10 1 

COS (8 

FIG. 7; Using the sampling method 45] for direct comparison 
between the experimental data and the expansion developed 
in [2^ . The dashed lines correspond to the uncertainty (la) 
in the determination of the coefficient combination i??" + /. 



efficients, as well as pinning down the only remaining 
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undetermined coefficient, ^3^^, using the present experi- 
mental event sample, is now a highly feasible plan for the 
immediate future. Thus another advancement has been 
made towards a complete characterization of the ampli- 
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unpolarized, in conjunction with the sampling method, 
we anticipate that all amplitudes can be determined in- 
dividually. For the future development of ChPT, detailed 
realistic constraints will be supplied. 
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